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13  ABSTRACT  I 

•''The  propagation  characteristics  for  three  model  magnetospheres  are  compared  for 
bydromagnetic  waves  in  the  micropulsations  frequency  range,  10~  rad/ sec 
1.0  rad/sec.  One  model  is  of  a  plane,  layered  medium  over  a  plane,  conduct ear on. 
a  second  model  is  one  composed  of  concentric  cylindrical  shells  of  plasma  surrounding  a 
conducting,  cylindrical  earth.  In  the  plane  model  the  magnetic  field  is  horizontal,  1  .t 
varies  with  height.  In  the  second,  the  stratic  mag  ietic  magnetic  field  is  parallel  to 
the  cylindrical  axis,  but  again  varies  with  height.  The  transmission  coefficient  for 
the  amplitude  of  the  magnetic  field  of  an  incident  plane  hydromagr.etic  wave  is  computed 
as  a  function  of  frequency  for  both  models,  and  the  results  are  compared.  The 
scattering  and  diffraction  effects  present  in  the  model  with  curved  surfaces  substanti¬ 
ally  reduce  the  transmission  coefficient  at  higher  frequencies.  Also,  differences  in 
geometric  shape  are  much  more  signnicant  than  any  reasonable  variation  of  other 
magnetospheric  parameters  within  either  model. 

A  third  model  utilizes  a  curved,  dipole- like  magnetic  field.  However, 
computations  are  not  complete  for  this  case. 


I 


TABLE  OF  CONTENTS 

PaRe 

I.  INTRODUCTION . 1 

II.  hAC-NETOSPHERlC  PROPAGATION  OF  HYPROMACNETIC  l,,.VES .  2 

A.  Introduction  .  .  . . 2 

B.  Specific  Magrietospheric  Models . 10 

C.  Results  of  Computations . 13 

III.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS .  17 

IV.  PERSONNEL . 18 

V.  REFERENCES . 19 

VI.  PUBLICATIONS . 20 


LIST  OF  FIGURES 


1.  The  variation  of  the  static  magnetic  field  B  -t  ebers/square 
meter)  and  mean  mass  of  positive  ions  (oxyger,  nlz-n  =  16) 

with  geometric  altitude  in  kilometers  ...  . . 21 

2.  Electron  number  density  K  (per  cubic  meter)  ,  s  a  function  of 

geometric  height  in  kilometers . . . ....22 

3.  The  magnetic  field  amplitude  transmission  coefficient  as  a 


function  of  angular  frequency  (radians/seccacl).  Solid  curve 

is  result  for  cylindrical  model;  broken  cu  «,  for  plane 

model;  collisions  have  been  neglected  in  T  >t'n  cases . 23 

4.  The  'ionospheric  magnetic  transmission  co"tficient’  t^  as 

calculated  by  Field  and  Greifinger  including  collisional 

effects  t.  "  as  calculated  in  this  paper.  Dashed  curve  is  the 

b 

ratio  t'/t,.  Abscissa  is  angular  frequency  in  radians/sccond.  24 

bo 

5.  The  magnetic  field  amplitude  transmission  coefficient  for  the 
two  cases:  P  represents  plane  model;  C,  cylindrical  model. 
Collisional  effects  included  in  both  cases.  Angular 

frequency  is  in  v.dians/second . ' . 25 


ABSTRACT 


The  propagation  characteristics  for  three  model  magnetospheres  are 

compared  for  bydroaagnetic  waves  in  the  aicropulsation  frequency  range, 

-4 

10  rad/sec  <  w  <  1.0  rad/sec.  One  model  5s  c-f  a  plane,  layered  median 
over  a  plane,  conducting  earth.  A  second  model  is  one  composed  of  concentric 
cylindrical  shells  of  plasma  surrounding  a  conducting,  cylindrical  earth. 

In  the  plane  model  the  magnetic  field  is  horizontal,  but  varies  with  height. 
In  the  second,  the  static  magnetic  field  is  parallel  to  the  cylindrical 
axis,  but  again  varies  with  height.  The  transmission  coefficient  for  the 
amplitude  of  the  magnetic  field  of  an  incident  plane  hydroE2g.:etic  wave  is 
computed  as  a  function  of  frequency  for  both  models,  and  thi  results  are 
compared.  The  scattering  and  diffraction  effects  present  in  the  model  with 
curved  surfaces  substantially  reduce  the  transmission  coefficient  at  higher 
frequencies.  Also,  differences  in  geometric  shape  are  much  more  significant 
than  any  resonable  variation  of  other  magnetospheric  parameters  within  either 
model. 

A  third  model  utilizes  a  curved,  dipole-like  magnetic  field.  However, 
computations  are  not  complete  for  this  case. 


I.  INTRODUCTiOK 


The  objective  of  the  research  supported  by  Air  Force  Contract 
F19628-67-C0152  vas  the  investigation  of  hydroaagnetic  phenomena  in 
the  magnetosphere  in  order  to  determine  the  possible  physical  mechanisms 
for  the  production  of  the  observed  geomagnetic  nicropulsations  vhose 
energies  are  confined  to  cooperatively  narrow  frequency  bandwidths. 
Essentially  all  the  work  performed  vas  directed  towards  understanding 
the  propagation  characterics  of  the  magnetosphere. 

This  work  vas  begun  January  2,  1967  and  concluded  vith  the 
termination  of  the  research  contract  on  December  31,  1970. 

The  work  described  in  this  report  is  primarily  the  theoretical 
determination  of  the  propagation  characteristics  of  low  frequency  hydro- 
magnetic  waves  in  three  model  magnetospheres :  (1)  a  plane,  layered 

medium  with  horizontal  magnetic  field;  (2)  a  cylindrical,  layered 
medium  with  static  magnetic  field  parallel  to  the  cylindrical  axis; 

(3)  a  cylindrical  model  with  static  field  normal  to  the  cylindrical  axis. 
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II.  MACSETOSPHERIC  PKOPAGAIIOX  OF  KYUSOlZAGSmC  WAVES 

A.  Introduction.  A  large  portion  of  the  work  covered  by  this 

report  was  published  in  a  paper  "A  Comparison  of  Jydronagnctic 
Wave  Propagation  in  Hagnetospheres  of  Plane  and  Cylindrical 
Geometries,"  Journal  of  Geophysical  Research,  Vol.  75,  pp.  3^1- 
3824,  1970.  Portions  of  the  nateri.'  1  in  sections  A,  B,  and  C 
of  this  chapter  are  quoted  from  this  paper 

The  mathematical  description  of  the  propagation  of 
hydroE3gnetic  waves  through  the  magnetosphere  is  complicated  by 
■any  factors:  (1)  the  mediu:..  is  inhomogeneous,  i‘s  properties 
varying  by  cany  orders  of  magnitude  over  distances  less  than  a 
wavelength;  (2)  the  static  magnetic  field  surrounding  the  earth 
causes  the  medium  to  be  anisotropic;  (3)  the  magnetospheric 
medium  is  dissipative;  (4)  the  parameters  tha*  characterize  the 
medium  are  not  known  experimentally,  except  approximately; 

X5)  the  geometry  of  the  magnetosphere  docs  not  lend  itself  to  a 
simple  mathematical  formulation  of  the  propagation  problem.  Because 
of  these  complications, a  realistic  magnetospheric  model  has  never 
been  used  in  propagation  calculations  at  hydronagnetic  wave  fre¬ 
quencies.  A  number  of  workers  have,  however,  reported  calculated 
result  based  on  approximate  models  [e.g.,  Field  and  Greifinger, 

1965,  1966;  Francis  and  Karpjus,  1960;  Greifinger  and  Greifinger, 
1965;  Jacobs  and  Watanabe,  1962,  Karplus  et  al.,  1962;  Prince  and 
Bostick,  1964;  Prince  et  al.,  1964]. 

A  typical  treatment  takes  into  account  some  of  the  factors 
listed  above.  A  'full  wave’  solution  to  the  propagation  equations 
can  be  affected  numerically,  even  for  the  inhomogeneous.,  anisotropic, 
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dissipative  medium,  if  tie  acgnetosphc-ric  parameters  are  functions 
of  only  one  Cartesian  coordinate.  In  such  an  analysis  the  earth 
becoees  an  infinite  half-space  with  a  plane  boundary,  the 
magnetosphere  is  layered  slab,  jor.d  the  fields  everywhere  are 
superpositions  of  plane  waves  E y  a3 loving  the  parameters  of  the 

medium  to  vary  with  height  in  a  sufficient' y  siople  fashion,  so:  *. 
workers  have  obtained  analytic  expressions  for  transnission 
coefficients  [Field  and  Greificger,  1965,  1966;  Gre: finger  and 
Greifinger,  1965].  In  general,  ch'-  tendency  has  beer,  to  try  coping 
with  the  first  four  difficulties  listed  previously,  out  to  restrict 
the  geometry  oi  the  mode*.  Sore  work  has  beep  done  or.  resonant 
codes  in  a  spherically  syrer-.'tric  mgnetos'here  'Radcski  and 
Corovillano,  1966;  Carcvillano  el  al. ,  196  >).  Ho  use  seems  to 
have  been  made  of  these  insults  in  predicting  transmission  charac¬ 
teristics,  however.  In  general,  the  analysis  of  magnerospheric 
propagation  cf  hydronagnctlc  waves  has  bean  o  dimensional  and  has 
made  use  of  plane  waves  and  plane  bounding  surfaces.  A  aiagnetospheric 
model  that  includes  curved  surface;  vcuV  also  include  scattering 
and  diffraction  effects  absent  in  th'  plane  layered  modejr. 

Thus  it  is  reasonable  to  attempt  a  solution  to  the  propagation 
problem  for  a  curved  (though  highly  artificial)  model  in  order 
to  evaluate  the  effects  of  curvature  and  hence  the  suitability  or 
accuracy  of  the  plane  laminar  models. 

Three  calculations  have  been,  attempted;'  two  have  been  com¬ 
pleted.  The  three  cases  are: 

1.  A  calculation  leading  to  a  soluti,".t  to  the  problem  of  the 
propagation  of  a  hydro:  ngnetic  disturbance  through  a  cyiindi  le.r  1  !y 
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symmetric  magnetosphere  of  a  plant  wave  incident  upon  it.  The 
continuously  variable  nature  of  the  medium  is  approximated 
by  a  number  of  homogeneous  cylindrical  shells  of  plasma.  The 
'dipole*  field  of  the  earth  in  this  mode]  varies  only  with 
distance  from  the  cylindrical  axis  and  is  parallel  to  the  axis. 

A  plane  normal  to  the  axis  of  this  model  is  thus  comparable  to 
the  equatorial  plane. 

2.  A  second  calculation  is  based  on  a  plane  model  of  the  magneto¬ 
sphere.  In  this  simpler  case,  the  sane  parameter  profiles  are 
used  as  in  the  cylindrical  model,  and  plane  layers  have  the 

same  thickness  and  height  as  the  corresponding  cylindrical  shells. 
When  the  transmission  coefficients  in  the  two  cases  are  compared, 
the  effects  of  curvature  of  the  magnet ospheric  surfaces  are 
revealed. 

The  exact  results  are  obtained  by  assuming  zero  collision 
frequencies  throughout  the  magnetosphere  in  both  cases.  However, 
the  effect  of  molecular  collisions  on  the  propagation  is  accounted 
for  approximately  by  making  use  of  the  ionospheric  transmission 
coefficient  computed  by  Field  and  G^eifinger  [1966]  to  modify 
thi  exact  no-co] li£ ion  results.  Although  this  composite  treatment 
is  not  a  rigorous  solution  to  the  propagation  equations  with 
collisions,  it  is  thought  to  be  a  reasonable  approximation  to  the 
proble.i.  In  addition,  this  approximate  solution  requires  con¬ 
siderably  less  conputirg  tine. 

3.  A  third  propagition  calculation  was  attempted  on  a  model 
magnetosphere  which  again  was  cyl indr i cally  symmetric.  However, 
in  this  third  case  the  magnetic  fi^ld  war.  chosen  tc  lie  in  the 
plane  normal  to  the  cylinder  a.  is  and  to  exhibit  a  "dipole"  . 
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characteristic.  Thus  the  direction  and  magnitude  of  the  static 
field  were  functions  of  position,  varying,  in  fact,  with  both  r 
and  the  radius  and  azimuthal  angle  in  cylindrical  coordinates 
This  considerably  more  complicated  magnetic  field  structure 
leads  of  course  to  a  much  more  difficult  mathematical  problem. 

At  the  date  of  termination  of  this  contract  the  solution  to 
the  propagation  problem  for  this  third  case  was  incomplete. 

A  description  of  the  progress  made,  however,  will  be  found  in 
section  C  beloi*. 

The  analysis  of  the  propagation  problem  for  the  first  two 
cases  proceeds  as  follows. 

Consider  a  homogeneous  plasma  medium,  corresponding  to  any 
one  of  th.  layers  in  the  magnetospheric  model.  Here,  Maxwell's 
equations  (nlcs  units) take  the  form 

7*E  =  p /cq  7*E  =  -jigII9/3t 

(1) 

7*H  «-  0  VxH  =  eg(3E/9t)  +  J 

For  time  dependence  of  the  form  e*ut:,  these  equations  become 

7*E  =  p/eg  V*e  =  “ipouH 

(2) 

*  0  7xH  =  icoegE  +  J 
c  itOCgK'E 

where  K  is  the  dielectric  coefficient  tensor.  If  the  plasma  is 
characterized  by  a  uniform  static  magnetic  field  %  parallel  to 
the  z  axis  of  a  Cartesian  system,  K  can  be.  written  as 


1 


K; 


-K. 


K 


0 

0 

K, 


(3) 


The  components  are  formed  from  the  basic  parameters  of  the 
medium,  X,  Y,  and  Z. 


X  *=  N  e2/ eom  Y  =  jeB0/m  ml 

s  s  s  s  1  s  ' 


Z 

s 


(4) 


where 

N  is  the  density  of  singly  charged  components  of  the 
s 

plasma. 

e  is  the  charge  cn  the  electron. 

is  the  ionic  mass  of  the  s*  component  of  the  plasma. 

v  is  the  collision  frequency  of  the  sth  component. 

s 

The  magnetospheric  plasma  will  be  treated  as  if,  at  any 
location,  there  existed  only  electrons  and  one  kind  of  positive 
ion,  which  has  a  mass  equal  to  the  average  mass  of  the  positive 
ions  at  that  location.  The  components  of  K  now  become 


K 


1  - 


X  8 
e  e 


(5) 


K 


x 


1Y  X 


e  e 


K 


iY.X. 
i  i 

ei^7 


(6) 

(7) 
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The  subscripts  i  and  e  refer  to  ions  and  electrons,  respectively. 

If  the  medium  is  homogeneous,  so  that  KA,  Kx,  and  Ky  are  independent 
of  coordinates,  and  if  the  fields  are  independent  of  z  (i.e., 
the  wave  propagates  normal  to  z) ,  Maxwell's  equations  can  be 
combined  to  form  the  wave  equations 

Vt2Ez  +  k02K/(E2  =  0  =  ?t2Ez  +  k22Ez  (8) 

^  Kj  +  K  2  _> 

Vt2Hz  +  k02  - — —  Hz  =  0  =  Vt21Iz  +  k!2Hz  (9) 

where  V^2  is  the  Laplacian  operator  in  the  pair  of  coordinates 
transverse  to  Bg  and  kg  is  equal  to  ui/c. 

Cylindrical  model.  The  solution  to  (9)  in  cylindrical 
coordinates  is 

00 

H  (r,<J>)  »  l  [A  J  (kjr)  +  B  Y  (kxr)]ein<f  (10) 

2  it  ii  (in 

na 

with’  J  ,  Y  Bessel  functions  of  the  first  and  second  kind  and 

n  n 

with  A  ,  B  arbitrary  constants.  If  A  and  B  can  be  evaluated 
n  n  n  n 

E  and  E.  can  be  determined  from  Maxwell's  equation.  Tnese 

r  $ 

three  field  components  form  the  'fast'  mode  which  can  propagate 

* 

normal  to  z.  The  other  possible  mode  conies  from  the  solution  to  (8) 

00 

E  (r  ,<J>)  -  l  [C  J  (kir)  +  DY  (k,r)]ein*  (11) 
n-  -OT 

The  two  modes  are  uncoupled,  since  propagation  is  normal  to  Bg. 

To  describe  the  fields  in  a  medium  that  consists  of  concentric 
cylindrical  layers  (each  layer  being  homogeneous  but  with  an  axial 

*  9 
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-u 

i 


static  magnetic  field),  we  begin  with  (10)  foi  one  mode  and 
derive  from  (2)  the  xesult 


,  rK,  3H 

F  =  1  _  _ z 

r  iweoCV  +  K  *)  r  3$ 


_ ±  _  If  _ ± 

3$  x  3r 


i  T  3H 

v  = _ 1 _  I  v  z  , 

*  io)£0 +  KxO  3r 


Sc  <] 

r  3*  j 


The  boundary  conditions  on  each  cylindrical  boundary,  r  -  a,  are 


V '  a~>  ’  Vr  -  “  > 


H  (r  -  a~)  =  H  (r  =  a  ) 
z  z 


Because  of  the  form  of  the  solution,  these  boundary  equations 
hold  independently  for  each  n  value  and  have  the  form,  for  example 


(kima)  +  BTim(k1ma)  =  Anm+1Jn(k1mfla)  +  Bnm+1Yn (k™+1a)  (16) 


n  n 


where  m  is  an  index  specifying  the  layer  in  which  the  field  is 
evaluated. 

In  the  present  problem,  one  boundary  will  be  at  the  magneto¬ 
pause.  Exterior  to  this  boundary  the  field  consists  of  an  incident 
plane  wave  and  a  scattered  cyl'.idrica]  wave.  The  form  of  the  H 
field  in  the  interplanetary  region  is  then,  for  the  fast  mode, 


H  (r,$)  =  H0  l  ineinSn(k!r)  +  £  F^^J^kjr)  (17) 

n=  -®  n=  -» 


■+  1Yn(kir)] 
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The  inneraost  boundary  to  the  asagaetospheric  region  is  at  the 
surface  of  the  earth.  The  earth  itself  is  assumed  to  be  a 
perfect  conductor,  so  that  the  electric  field  vanishes  for 
r  <  Kg.  The  boundary  condition  is  then 

Vv** =  0  (18) 

The  conditions  above  (equations  16  and  18)  provide  just  the 
peeper  number  of  conditions  to  allow  a  solution  to  be  formed 
for  each  undetermined  constant  in  the  field  expression  (equations 
10,  12,  and  13,  for  example).  In  a  model  of  N  layers  (counting 
both  the  earth  and  the  interplanetary  medium),  the  number  of 
unknowns  is  2N  -  3.  However ,  because  these  unknowns  are  complex 
there  are  4N  ~  6  values  to  be  determined,  and  hence  4N  -  6  real 
equations.  Solution  to  the  set  of  4N  -  6  equations  determines 
the  coefficients  in  the  field  expansions  for  one  n  value  and  one 
frequency. 

Plane  model.  For  propagation. normal  to  the  plane  boundaries, 
equation  9  has  the  solution 


H2(x>  -  4-  (19> 

where  Hzj  and  hZ2  are  the  complex  amplitudes  of  the  plane  •  vis 
propagating  in  the  positive  and  the  negative  x  directions,  respectively. 
The  corresponding  electric  field  is 


-k,K,H  , 

E  (x)  - - tH- 

y v  weo(K  ' 


+  K2) 

X  7 


iklX 


toco  (K,  1  +  K/) 


-ikix 


(20) 


Ex(x)  =  (Kx/Ka)Ky(x) 


(21) 

t 
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The  boundary  cond  \i<«ns  on  the  plane  boundaries  (x  is  constant) 
are  that  and  be  continuous.  Again,  4S  -  6  equations  are 
solved  for  the  AS  -  6  unknown  amplitude  values.  Only  the  fast 


mode  solution  is  considered. 


Specific  Magnetospbcric  Models 

1.  Cylindrical  Magnetosphere.  The  earth-*sagnetospl»ere-interplanetary- 
space  system  is  represented  in  vhc  curved  oodel  by  a  set  of  twenty- 
two  concentric  cylinders.  The  innemost,  ».hich  is  the  earth,  is 
a  perfect  conductor.  The  next  in  order  of  increasing  radius  is 
the  neutral  nonconducting  atmosphere.  The  following  nineteen 
cylindrical  layers  represent  the  ionosphere  magnetosphere.  The 
last  'layer*  is  bounded  on  the  inside  by  the  magnetopause,  but  it 
has  an  infinite  outer  radius.  Each  layer  is  considered  to  be 
homogeneous,  with  parameters  as  indicated  in  Table  1. 

The  values  assigned  to  the  parameters  for  each  layer,  together 
with  the  radii  appropriate  to  each  layer,  were  decided  as  follows. 
First,  a  table  of  mean  molecular  weight,  electron  number  density, 
and  geomagnetic  field  intensity  was  compiled  for  fifty-five 
geometric  altitudes  between  60  and  61,000  km.  These  valuer  were 
taken  from  the  literature  [Prince  et  al.,  1964;  Matsushita  and 
Campbell,  1967].  Smooth  curves  of  the  three  parameters  are  shown 
in  Figures  1  and  2.  In  this  model  no  consideration  was  given  to 
collisions  of  ions  with  neutrals  or  with  other  ions.  The  consequences 
of  this  neglect  will  be  discussed  later.  Values  were  chosen  to 
represent  approximately  sunspot  minimum,  daytime. 

The  values  of  the  log  of  the  Alfven  velocity  V  calculated  at 

A 

the  fifty-five  data  points  were  plotted  against  log  R  and  a  smooth 
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curve  teas  drawn  through  the  points.  The  ionosphere-magnetosphere 
was  divided  into  nineteen  regions;  their  boundaries  were  chosen 
simply  by  inspection,  hit  in  such  a  vay  as  to  fit  the  cane 
ns  closely  as  possible,  as  judged  by  ere.  The  set  of  radii  thus 
chosen  were  inposed  on  smooth  curves  of  Kj^  and  values  for 
each  region  selected.  The  mother  of  layers  cboseu  (twenty- two) 
was  a  cooprosise  between  an  accurate  fit  and  a  reasonable  anount 
of  computer  tine. 

2.  Plane  Magnetosphere.  The  sase  parameters  were  used  to  describe 
the  plane  layered  node.  Each  radius  listed  in  Table  1  now  becoees 
a  height  x  above  the  center  of  the  plane  conducting  slab  that 
represents  the  earth. 


« 
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TABLE  1-  Piranttfrs  of  the  Tw«ty-TV?-Lircr  Kodel 


Layer 

Badii, 

meters 

Tj,  if /sec 

** 

(a*  =  0.1/sec) 

K 

(u  *  0.1/sec) 

1 

0 

6.37 

X 

10* 

•  •• 

•  •  • 

... 

2 

6.37 

X 

10* 

6.43 

X 

10* 

3.00  x  10* 

0.00 

l.ce 

3 

6.43 

X 

10* 

6.47 

X 

10* 

2.00  x  10* 

4.00  x 

10 

5-00  x  10* 

4 

6.47 

X 

10* 

6.51 

X 

10* 

4.40  x  10* 

5.40  x 

10* 

6.00  x  10s 

5 

6.51 

X 

10* 

6.61 

X 

10* 

2.90  x  105 

1.10  x 

103 

1.40  x  10s 

6 

6.61 

X 

10* 

6.83 

X 

10* 

2.30  x  105 

1.40  x 

103 

2.15  x  10* 

7 

6.83 

X 

10* 

7.17 

X 

10* 

4.00  x  105 

5.00  x 

107 

9.00  x  105 

8 

7.17 

X 

10* 

7.67 

X 

10* 

1.00  x  10* 

6.00  x 

10 

1.10  x  10s 

9 

7.67 

X 

10* 

8.27 

X 

10* 

2.30  x  10* 

3.00 

1.70  x  10* 

10 

8.27 

X 

10* 

1.08 

X 

10* 

4.80  x  10* 

4.60  x 

10-1 

4.80  x  103 

11 

1.08 

X 

107 

1.18 

X 

107 

3.85  x  10* 

1.25 

7.5C  x  10* 

12 

1.18 

X 

107 

1.34 

X 

107 

3.25  x  10* 

2.30 

1.00  x  10* 

13 

1.34 

X 

107 

1.52 

X 

107 

2.70  x  10* 

5.30 

1.60  x  10“ 

14 

1.52 

X 

107 

i.88 

X 

107 

2.05  x  10* 

1.70  x 

10 

2.60  x  10* 

13 

1.88 

X 

n7 

2.24 

X 

:o7 

1.52  x  10* 

6.00  x 

10 

4.60  x  10* 

16 

2.24 

X 

107 

2.54 

X 

107 

1.70  x  10* 

3.80  x 

10 

5-80  x  10* 

17 

2.54 

X 

107 

2.96 

X 

107 

3.50  x  10* 

2.20  x 

10 

5.40  x  70 

18 

2.96 

X 

107 

3.44 

X 

107 

7.80  x  10* 

3.30  x 

10 

4.50  x  10* 

19 

3.44 

X 

107 

4.64 

X 

107 

2.85  a  10* 

1.30  x 

10‘ 

3.40  x  10* 

20 

4.64 

X 

107 

6.20 

X 

107 

1.50  x  10* 

9.20  x 

102 

4.50  x  10* 

21 

6.20 

X 

107 

6.80 

X 

107 

1.45  x .10* 

1.25  x 

103 

5.00  x  10* 

22 

6.80 

X 

107 

m 

2.12  x  105 

4.31  x 

10* 

2.58  x  10* 

Layer  1  repre  'nts  the  earth;  layer  2,  neutral  atmosphere;  layers  3-21,  the 
ionosphere-magnetosphere;  layer  22,  interplanetary  space. 
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1.  Calculations  of  fcgcetic  Field  Amplitude  Traasatijsico 
Coefficient  f or  First  Two  Cases.  Use  cylindrical  model  lecd-i 
to  a  set  of  eifhty-tuo  simultac^eus  real  linear  algebraic 
•^uatiens  for  the  coefficients  A  *  1  *,  and  F  in  the  field 
expressions  (ecaations  10  and  17);  tite  superscript  designates 
the  nth  layer.  One  value  of  u  was  chosen,  and.  the  set  of 
equations  was  solved  for  each  integer  n  in  the  interval  -  10  < 

n  <  +  10.  The  values  for  A  2  and  B  2  were  then  substituted  into 

*=  n  n 

(10),  and  the  magnitude  ^  iron  2  to  100.  The  trend  of  the 
ratio  to  increase  with  increasing  frequency  is  apparent.  The 
resonant  frequencies  for  both  cases  fell  within  the  range  of 
continuous  sinusoidal  oscillations  observed  at  geostationary 
altitudes  (0.C2  to  0.12  radian/sec)  and  interpreted  as  nagneto- 
sphe.ic  resonance  standing  waves  [Cunnings  et  al..  1969]. 

The  low  frequency  linits  for  the  two  cases  differ  by  a  factor 
of  2.  This  result  could  have  been  anticipated,  since  at  these 
frequencies  the  magnetosphere  has  little  effect  on  the  wave. 

Thus,  as  u  -*  0,  the  comparison  is  between  the  reflection  at  an 
infinite  conducting  plane  and  the  scattering  by  a  conducting 
cylinder.  The  fields  calculated  in  these  cases  exhibit  a  factor 
of  2  in  fh  ■  ratio  of  the  amplitudes  of  the  magnetic  intensity  at 
the  surfaces. 

Field  and  Creifinger  {1966]  defined  their  ’ionospheric  magnetic 
transmission  corf f icicnt , '  t„,  as  the  ratio  of  the  total  magnetic 

D 

field  at  the  bare  of  the  ionosphere  to  the  field  of  the  dovneo:  ing 


u 
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fort  of  the  sigaal  at  the  tcj  of  tW  i»»jVre.  Tbe  top  of  the 
imo$pi»re  V2t  uba  to  k  430  b  zkv*  tie  surface  for  daytime, 
sanspct  naglacai  coeifticss;  tie  bar*  of  tk  josc^kre,  to  be 
S3  ke  hip.  Tbe  node 1  of  vildt  .ky  S2k  ese  has  a 

constant  Alves  speed  is  this  rtfisa,  ¥  *  2 GO  km/sec,  a  coast amt 

in  cydotroo  freyseacy  of  190  uiia^/sft,  and  ioa  zad  electros 
collision  frtqccacics  proport local  to  each  other  and  both  exponential 
functions  of  height.  Obviously  the  presect  model  is  mo:  identical 
to  that  of  l’ield  and  Creifinger,  eves  aside  front  collision  effects. 

Kever the less ,  a  comparison  of  their  ccajwled  with  as  ionospheric 
magnetic  transatissioa  coefficient  calcslatc-d  os  the  basis  of  the 
present  plane  layered  node 1  should  give  scute  indications  of  the 
attenuation  resulting  from  collisions-  Such  a  comparison  is 
Illustrated  in  Figure  4,  where  t  '  is  the  corresponding  transmission 
coefficient  derived  f ran  the  present  work.  The  evaluation  of  t  ' 
was  oade  by  cosparing  the  amplitude  of  the  nagnetic  field  at  the 
surface  with  the  anplitude  of  an  incoming  plane  wave  above  the 
boundary  at  6.83  x  10  neters.  or  460  fcs  above  the  surface.  Ionospheric 
reflections  of  the  incident  wave  causes  the  gradual  decrease  in  t  ' 
ns  the  frequency  is  increased,  until  the  layered  red i urn  becomes 
resonant  at  above  6  radians/sec.  The  t  curve  decreases  rapidly 
because  of  absorption,  and  sc  shows  no  resonance.  The  ratio  of  the 
two  curves  t /t  '  is  shown  in  the  broken  cur’-c  in  Figure  4. 

To  cialce  use  of  the  Field  and  Gre.' finger  result,  the  spectra 
of  the  plane  and  cylindrical  magneto spheres  verc  multiplied  by  the 
ratio  t  '/*  .  The  application  of  this  quantity,  based  on  ?  one— 
dimensional  calculation,  to  the  cylindrical,  model  sccr.s.  justified 

'*  # 
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dace  the  t’kicbcst  of  tk  uoo^kric  layer*  iarclrd  i*  small 
la  conpzrisMt  with  the  radius  of  temtare  of  the  booodiag 
nrfacts,  aanelj  4W  Li/637?  ia,  or  dov:  71.  The  final  spectra 
lacWicg  the  collision  effects  are  shown  is  Tip, re  5.  The 
ratio  of  these  two  corses,  f  for  the  plane  layered!  jodel  and 
C  for  the  cylindrical  one,  is  of  course  still  the  same  as  In 
figure  3. 

One  difference  in  the  two  traasaissioc  carves,  in  addition  to 
the  depressed  transmission  in  the  cylindrical  case  brought  about 
by  scattering,  is  the  shift  in  resonant  frequencies.  The  first 
resonance  comes  at  a  higher  frequency  in  the  cylindrical  model 
than  in  the  plane  node.  Examination  of  the  field  amplitude  through¬ 
out  the  medium  indicates  that  this  effect  is  a  result  of  multiple 
internal  reflections  in  the  plane  model,  which  lead  to  a  longer 
effective  path  length  for  the  wave  propagating  normal  to  the 
plane  bounding  surfaces.  The  higher  resonant  frequencies  are 
harmonically  related  tc  the  first  resonance  in  the  cylindrical 
Model,  but  their  amplitudes  are  rapidly  decreasing,  owing  to  the 
scattering  at  each  curved  boundary.  The  relationship  between 
resonant  frequencies  is  more  complicated  in  the  plane  model  because 
of  aultiplc  reflection  and  individual  layer  resonance  occurring 
at  frequencies  above  0.35  radian/sec. 

The  sensitivity  of  the  results  tc  the  details  of  the  models  was 
investigated  by  varying  the  nagnetospheric  para-.elcrs  >ver  a  rather 
wide  range  of  values.  In  general,  the  spectra  are  oily  slightly 

sensitive  to  such  changes.  To  illustrate  this  insensitivity,  a 

*  $ 
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calccUtjoc  vats  wade  for  both  wodcls  is  vbicb  taut  satire  najKtfr- 
Sfbe re  loowpfetre  was  replaced  by  *  single  hmogmoers,  isotropic 
layer.  The  dielectric  constant  of  this  layer  vas  adjusted  to 
give  the  see  propagation  time  as  in  the  detailed  sodcls.  The 
result  in  the  cylindrical  case  is  a  spectrum  displaying  resonances 
that  coincide  with  the  first  fr«  restores  in  the  nere  elaborate 
■odd  but  haring  a  transmission  coefficient  that  falls  off  more 
slowly  with  frequency.  The  fact  that  the  transmission  coefficient 
is  not  strongly  coupled  to  the  specific  details  of  the  aajnetospheric 
■odd  has  tvo  implications:  First,  the  approximations  used  in  the 
daculatioes  for  the  cylindrical  model  tend  to  be  justified;  second, 
the  transmission  coefficient  for  the  magnetic  field  amplitude  is 
seen  to  be  ouch  core  strongly  dependent  on  the  geometry  of  the-  magnet  o- 
Sjheric  model  than  on,  for  example,  the  details  of  parameter  profile*  . 
2.  Progress  in  Tvo-diner.sional  Dipole  Case.  A  "two-diaensieual 
dipole"  field  vas  assumed: 


pH 


^(r»^)  =  ~7  sin  f  r  +  c°s  (•  $) 


ir 

The  wave  equations  for  the  wave  field  components  are  strongly 
coupled  by  this  static  field  structure-  A  nethod  of  decoupling 
lccally  vas  found,  which  in  principal  would  allow  solution  to  the 
wave  equations  in  a  snail  volune-  However,  the  boundaries  for 
these  snail  volumes  were  such  that  the  boundary  condition  could 
not  easily  be  ratified.  A  corprorise  static  field  structure  was 
next  employed.  This  ar.bient  field  was  radical  at  high  "lat ' tudes" 
(^r  <_  p  <  — and  was  azimuthal  at  low  latitudes  (-  ^  <  p  <  *■  ). 
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This  static  field  has  wsit  chancteristics  of  the  dipole  'ield, 
but  again  bctadzr)  cocdftiois  could  not  be  solved.  Finally,  it 
was  found  that  in  a  thin  cylindrical  shell  (of  thickness  t  where 
t  «  r,  and  r  =  inner  radius  of  the  shell)  the  wave  equation 
could  be  decoupled.  For  thi;-  case  the  wave  equation  and  the 
dispersion  relation  can  be  solved  numerical ly,  and  the  boundary 
conditions  satisfied.  However ,  numerical  calculations  aro  not 
complete.  Thus,  the  evaluation  of  the  propagation  of  h-u  waves 
through  a  "two  diinsioaal  dipole"  field  in  a  cylindrical  cagneto- 
sphere  was  weakened  to  the  determination  cf  propagation  through  a 
thin  cylindrical  shell  of  plassa  with  an  azicuthal  field.  This 
problen,  while  soluble,  was  not  completed  in  the  contracting 
period. 


m.  sukmaky.  conclusions,  jgp  rkco^:f..'datio:;s 

A.  Suraary  and  Conclusions.  A  conparison  of  the  magnetic  field 

amplitude  transmission  coefficient  for  two  aodel  magnetospheres, 
one  cf  cylindrical  geometry  and  the  other  plane  geometry,  reveal, 
a  significant  difference.  The  curved  magnetosphere  scatters 
and  diffracts  an  incident  wave,  and  so  transnits  less  energy  to 
the  surface.  The  ratio  of  the  transmission  coefficients  in  the 
two  cases  varies  fron  2  to  about  100  ovei  the  hydromagnatic  range 
of  frequencies.  The  lowest  resonant  frequencies  are  seen  to  depend 
largely  on  the  size  of  the  magnetosphere  and  average  velocity  of 
hydromagnetic  waves.  However,  a  plane  layered  model  may  produce 
more  complicated  resonant  behavio.  because  of  multiple  internal 
reflections.  A  fault  of  a  onc-dii-cnsion-il  treatment  is.  its  over¬ 
emphasis  of  the  effect  of  the  earth  on  magnetospheric  propagation. 
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a  fault  iabereat  in  the  place  node.  Previous  calculations  that  were 
based  oo  a  codel  involving  plane  layers  and  that  indicated  ampli¬ 
fication  effects  at  discrete  frequencies  should  be  reexaain ed  if 
comparison  between  interplanetary  field  and  surface  field  amplitudes 
is  desired. 

The  evaluation  of  propagation  characteristic  in  a  magnetosphere 
with  curved  static  field  vas  attempted,  but  not  completed.  Such  a 
calculation  would  suggest  possible  lati  udinal  variation  in  surface 
wave  field  intensities  due  to  the  guiding  tendency  of  the  ambient 
field. 

B.  Recogrendations.  It  is  rec err ended  that  additional  effort  be  put 
into  the  curved  field  calculations.  Scattering,  diffracting,  and 
guiding  effects  in  this  case  should  m2kc  the  propagation  characteristics 
differ  considerably  froa  the  two  previous  calculations,  and  should 
provide  an  indication  of  the  latitutional  variations  in  surface 
field  intensities  to  be  expected  from  h-n  waves  propagated 
through  the  magnetosphere. 

IV.  PERSONNEL 

The  personnel  employed  under  this  contract  werr  Dr.  A.  W.  Jenkins,  Jr., 

Mr.  Hector  Dinas,  and  Mr.  William  Lee.  Mr.  Dinas  and  Mr.  Lee  were  employed 
as  graduate  research  assisr-nts.  Dr.  Jenkins  was  principal  investigator. 


19 


V.  REFERENCES 


1.  Allis,  W.  P. ,  S.  J.  Buchsbaua,  and  A.  Beers,  Waves  in  Anisotropic 
Plasmas,  M. I.T.  Press,  1963. 

2.  Carovillano,  R.  L. ,  H.  R.  Radoski ,  and  J.  F.  McClay,  Foloioal 

hydroEagaetic  plasmaspheric  resonances,  Phys.  Fluids,  9 ,  1860- 
1864,  1966.  *  "" 

3.  Cumeiegs,  W.  D. ,  R.  J.  O'Sullivan,  and  P.  J.  Coleman,  Jr.,  Standing 
Alfven  waves  in  the  ionosphere,  J.  Geophys.  Res.,  74,  778-793,  1969. 

4.  Dimas,  H.,  "Alfven  Wave  Propagation  in  a  Cylindrical  Magnetosphere”, 
Thesis,  l.orth  Carolina  State  University,  1969. 

5.  Field,  E.  C. ,  and  C.  Greifinger,  Transmission  of  geomagnetic  micro- 
pulsations  through  the  ic nosphere  and  lower  exosphere,  J.  Geophys. 

Res.,  70,  4885-4S99,  1965. 

6.  Field,  E.  C. ,  and  C.  Greifinger,  Equatorial  transmission  of  geomagnetic 
■icropulsations  through  the  ionosphere  and  lower  exosphere,  J.  Geophys. 
Res. ,71,  3223-3232,  1966. 

7.  Francis,  W.  E. ,  and  R.  Karplus,  Hydroniagr.atic  waves  in  the  ionosphere, 

J.  Geophys.  Res.,  65,  3593-3600,  1960. 

8.  Greifinger,  C. ,  and  P.  Greifinger,  Transmission  of  micropulsations 
through  the  lower  ionosphere,  J,  Geophys.  Res.,  70,  2217-2232,  1965. 

9.  Jacobs,  J.  A.,  and  T.  Watanabe,  Propagation  of  hydromagnet ic  waves 
•  in  the  lower  exosphere  and  the  origin  of  short  period  geomagnetic 

pulsations,  J.  Atmos.  Terr.  Phys.,  24,  413-434,  1962. 

10.  Jenkins,  A.  W.  Jr.,  "A  Comparison  of  Hydromagnetic  Wave  Propagation  in 
Magnetospheres  of  Plane  and  Cylindrical  Geometries,”  J.  Geophys.  Res., 

75,  3817-3824,  1970. 

11.  Karplus,  R. ,  W.  E.  Francis,  and  A,  J.  Dragt,  The  attenuation  of 
hydromaguetic  waves  in  the  ionosphere,  Planet.  Space  Sci. ,  j),  771-785, 

1962. 

12.  Matsushita,  S. ,  and  W.  H.  Campbell  (Ed.),  Physics  of  Geomagnetic 
Phenomena,  2  vol. ,  Academic  Press,  New  York,  1967. 

13.  Prince,  C.  E. ,  F.  X.  Bostick,  Ionospheric  transmission  of  transversely- 
propagated  planewaves  at  micropulsation  frequencies  and  theoretical 
power  spectrums,  J.  Geophys.  Res.  ,  69,  '’213-1234,  1964. 

14.  Prince,  C.  E. ,  F.  X.  Bostick,  and  II.  W.  Smith,  A  study  of  the  transmission 
of  plane  hydromagnetic  waves  through  the  upper  atmosphere,  Univ.  Texas 
Elec.  Eng.  Res.  Lab. ,  Rep. _ 134,  July  1964. 

15.  Radoski,  H.  R. ,  and  R.  L.  Carovillano,  Axisyminetric  plasmaspherc 
resonances.:  Toroidal  mode,  Phys.  Fluids,  9,  285,  1966. 


20 


VI.  PUBLICATION'S 

Two  Publications  resulted  from  this  contract: 

1.  A.  W.  Jenkins,  "A  Comparison  of  Hydromagnetic  Wave  Propagation 
in  Magnetospheres  of  Pland  and  Cylindrical  Geometries",  J. 

Geophys.  Res.,  75,  3817-382 4,  1970 

2.  Hector  Dimas,  "Propagation  of  Hydromagnetic  Waves  in  a  Cylindrical 
Magnetosphere,"  Thesis,  N.  C.  State  University,  1969. 


